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Abstract 

Periodic crustal deformation associated with repeated strike 
slip earthquakes is computed for the following model: A depth L 
(<fl) extending downward from the Earth's surface at a transform 
boundary between uniform elastic lithospheric plates of thickness H 
is locked between earthquakes. It slips an amount consistent with 
remote plate velocity V pi after each lapse of earthquake cycle time 
T cy . Lower portions of the fault zone at the boundary slip 
continuously so as to maintain constant resistive shear stress. 
The plates are coupled at their base to a Maxwellian viscoelastic 
asthenosphere through which steady deep-seated mantle motions, 
compatible with plate velocity, are transmitted to the surface 
plates . The coupling is described approximately through a 
generalized Elsasser model. We argue that the model gives a more 
realistic physical description of tectonic loading, including the 
time dependence of deep slip and crustal stress build-up throughout 
the earthquake cycle, than do simpler kinematic models in which 
loading is represented as imposed uniform dislocation slip on the 
fault below the locked zone. Parameters of the model are chosen in 
accordance with seismic and geologic constraints and to fit the 
apparent time-dependence, throughout the earthquake cycle, of 
surface strain rates along presently locked traces of the 1857 and 
1906 San Andreas ruptures. We find that prediction based on the 
resulting parameters compare reasonably to data on variations of 
contemporary surface strain and displacement rates as a function of 
distance from the 1857 and 1906 rupture traces, although the data 
is generally affected by asymmetry relative to the fault and by 
adjacent fault strands. Specifically, we fix V. = 35 mm/yr, T = 
160 yr and L = 9 to 11 km as a representative earthquake nucleation 
depth with a 2 km allowance for possible upward motion of the 
locked zone border during the earthquake cycle . We then find that 
the geodetic data is described reasonably, within the context of a 
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model that is locally uniform along strike and symmetric about a 
single San Andreas fault strand, by lithosphere thickness H = 17 to 
25 km and Elsasser relaxation time t r = 10 to 16 yr. We conclude 
that the asthenosphere appropriate to describe crustal deformation 
on the earthquake cycle time scale lies in the lower crust and 
perhaps crust-mantle transition zone, and has an effective 
viscosity between about 2 x 10 18 and 10 19 Pa-s, depending on the 
thickness assigned to the asthenospheric layer. 
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Introduction 

The northern part of the San Andreas fault in California was 
last broken by the great 1906 earthquake rupturing a 450 km 
segment. The southern part was broken last in 1857 with a rupture 
length of 350 km. Both of these fault segments have been locked 
since their last earthquake, with attendant strain accumulation 
over time. These locked segments are characterized by low levels 
of seismicity (see e.g., Carlson et al, 1979). In contrast, the 
San Andreas fault in central California appears to be creeping so 
that no long term strain accumulations occurs. 

The nature of strain accumulation at a strike slip plate 
boundary has been discussed by various authors. In general, the 
strain rates are elevated near the fault and this is presumed to be 
due to deep aseismic slip or distributed shear flow along the 
downward continuation of the fault plane. The near surface fault 
locking and deep aseismic deformation is thought to be a result of 
material property contrasts. Near surface, the fault may generally 
be characterized as elastic and brittle. At greater depth, the 
material may be undergoing plastic shear flow due to the high 
temperature and pressure. Thus deformation would be expected to be 
aseismic and to accummulate continuously at depth while the upper 
crust has to accommodate relative movements between the North 
American Plate and the Pacific Plate by seismic faulting. Sibson 
(1982) and Meissner and Strehlau (1982) gave an account of the 
changing material properties with depth and explained the shallow 
confinement of seismicity on the San Andreas Fault in terms of a 
transition from brittle friction to ductile creep. Tse and Rice 
(1986) showed in an elaboration of this concept that the 
temperature and hence depth variation of fault slip constitutive 
properties, when incorporated into a mechanical model for plates 
joined at a transform boundary, led to predicted depth variable 
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slip consistent with a shallow effectively locked zone that 
undergoes great seismic slips and with continuous slip accumulation 
below. Over the last several decades, strain changes have been 
monitored by geodetic networks at several locations along the 
fault. These data provide a basis for better understanding of the 
physical processes governing the behavior of the fault and, 
together with seismic and geologic data, are a source of possible 
constraints on locked zone depth, lithosophere thickness and 
rheological parameters . 

Our goal here is to develop a simple physical model of the 
earthquake stressing process along the locked segments of the San 
Andreas fault that is compatible with such data. While we do not 
discuss here events leading to great earthquake instability, we 
expect that the improved description of the fault region geometry 
and rheological properties will be useful to more realistically 
representing the earthquake loading process in necessarily more 
complicated crustal scale instability models. Such models have 
been based on shear crack concepts (Dmowska and Li, 1982; Li and 
Rice, 1983; Tse et al., 1985; Li and Fares, 1986), fault surface 
slip-weakening response (Stuart and Mavko, 1979; Stuart, 1984/85; 
Stuart et al, 1985) or its generalization in the manner of slip 
rate and slip history dependent friction (Mavko, 1980; Tse and 
Rice, 1986) . 

Savage and Burford (1973) suggested the modelling of 
interseimic surface strain rate profiles near strike slip faults by 
means of a buried screw dislocation in an elastic half-space. The 
technique has been employed widely (e.g., Prescott et al., 1979; 
McGarr et al., 1982; Savage, 1983; King and Savage, 1984). In that 
approach there is assumed to be no slip at the transform margin 
within a presently locked seismogenic depth range, analogous to 
what we denote as L here, whereas a spatially uniform slip rate at 
the relative plate velocity V pl is imposed at the margin everywhere 
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along the downward continuation of the locked zone. Savage and 
Bur ford (1973) recognized this imposed motion on the plate boundary 
as a convenient simplification of a more complicated driving 
mechanism in which "the plates slip past one another in response to 
shear stresses that presumably orginate from the drag of mantle 
convection currents upon the bottom of lithospheric plates". 
Evidently, the motion at the transform margin below the locked zone 
would be a calculated response, rather than as imposed condition, 
in a more complete model of driving by deep-seated mantle motion 
and its coupling to the surface plates. While the modelling which 
we develop here is still severely simplified, it does nevertheless 
come closer to embodying the driving mechanism just described. 

Turcotte and Spence (1974) analyzed the near-fault surface 
deformation by means of a two dimensional elastic edge crack model, 
like that in Fig. lb. The crack faces represent the sliding 
portion of the plate margin and are assumed to slide under 
resistive shear stresses below the locked zone which remain uniform 
in time, with loading from the remote edges of the strip. Since 
the interest is in changes of displacement under ongoing loading, 
this constant stress condition is equivalent to treating the crack 
surfaces as freely slipping. Their analysis provided an 
appropriate assessment of the effect of the depth of locked zone on 
the profile of surface strain or displacement, particularly on the 
ratio between the near fault value to that at a point remote from 
the fault (where stress and strain alterations are presumably 
uniform through the thickness) . The Turcotte-Spence crack model 
for the deeper fault zone, with its constant stress boundary 
condition, provides an alternative and presumably more realistic 
treatment of that zone than imposing a spatially uniform 
dislocation accumulating at a constant rate. That is, if the 
deeper fault zone deforms viscously with a strongly non-linear 
stress dependence as expected, then the local shear stresses will 
vary only modestly over appreciable changes in slip rate and are 
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sensibly approximated as being uniform in time except in the 
co-seismic and short-term postseismic time intervals. Also, the 
resulting crack surface slip distribution tapers gradually to zero 
at the base of the locked zone, eliminating the unrealistic slip 
discontinuity of the uniform dislocation model. 

Despite the attractive features of the crack model, its level 
of development by Turcotte and Spence remains of limited 
application because they assumed for simplicity in their modelling 
that the surface plates were decoupled from at least the nearby 
mantle below. Thus, they could only load the system by remotely 
applied forces and could not directly relate the loading to ongoing 
plate motion. 

We remove these deficiencies so as to account approximately 
for the coupling of the lithosphere to a viscoelastic asthenosphere 
and to allow the overall lithosphere-asthenosphere system to be 
loaded by ongoing deep seated mantle motions that are compatible 
with remote plate velocities. As will be seen, this more complete 
analysis of the edge cracked model allows us to fit approximately 
the principal types of geodetic data now available, namely data on 
the time dependence of near fault deformation over the earthquake 
cycle and the variation of comtemporary deformation rates with 
distance from the fault trace. 

The above mentioned works by Savage and Burford (1973) and 
followers, and by Turcotte and Spence (1974) , being purely elastic, 
do not account for observed changes in strain rates over a complete 
earthquake cycle. Thatcher (1975) proposed a qualitative model for 
strain accumulation and release of the 1906 San Francisco 
Earthquake. Essentially, the fault was locked at the upper 10 km 
which was the focal depth. Tectonic plate movement drove the plate 
boundary deformation, which was localized along the fault due to 
aseismic sliding occurring below the locked zone. The accumulated 
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strain was released by the 1906 earthquake, rupturing upwards 
through the locked portion. Post-seismically, the deeper part of 
the fault slid more rapidly, being driven by the stress shed onto 
it and the asthenosphere below by the earthquake. Thus, the 
postseismic surface strain rate was higher, but decreased gradually 
with time. For example, the (engineering) strain rate near the 
fault was about 2.4 x 10“ 6 yr -1 for the thirty years following the 
earthquake but was only 0.6 x 10 -6 yr -1 since then (Thatcher, 1975). 
The time-dependence of crustal deformation near a fault may have 
two major sources: inelastic relaxation of the fault zone material 
below the seismogenic layer, and coupling between the elastic 
lithosphere and the viscoelastic asthenosphere. Thatcher (1983) 
modelled the aseismic deep slip by means of an elastic half space 
in which postseismic transient slip with exponential time decay is 
imposed kinematically as a spatially uniform dislocation. A 
Nur-Mavko model, in which an elastic layer is coupled to a 
viscoelastic half space is used to model the latter effect. This 
work provided important insight on the nonlinear nature of strain 
accumulation with time over an earthquake cycle. The strain rate 
data from northern and southern California (Thatcher, 1983) argues 
strongly in favor of this nonlinear accumulation. Our work here 
addresses the same data, but without the kinematic imposition of 
motion directly beneath the seismogenic zone. Rather, the deep 
fault walls move, in our modelling, in response to steady mantle 
motions, as transmitted through a viscoelastic asthenosphere, so as 
to maintain the constant resistive stress boundary condition. 
Turcotte et al (1984) addressed the time dependent straining with a 
multi-layer model, having an intracrustal as well as a deeper 
asthenosphere, that is somewhat in the spirit of the present 
approach, but did not include the strain concentrating effect of 
slip on deep fault walls within the elastic surface plate. 
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The Model 

We build-in the viscoelastic asthenospheric coupling between 
the surface plate and deeper mantle motions by use of the 2-D 
generalized Elsasser model (Rice, 1980; Lehner et al, 1981) . This 
results in a differential equation which is phrased in terms of 
thickness average stress and displacement. Further, to represent 
the edge cracked geometry with slipping crack walls at the plate 
margin. Fig. lb, we require that the local displacement at the 
margin be related to the net shear force transmitted across it 
exactly as in the Turcotte-Spence analysis . 

Thus, considering as in Fig. la, an elastic lithospheric plate 
underlain by a viscoelastic asthenosphere, for the strike slip 
environment of a long fault with uniform conditions along strike 
equilibrium requires that 

5a xy /3y = x x /H (1) 

Here cJ X y(y,t) is the thickness average of the fault-parallel 
in-plane shear stress over plate thickness H and T x (y,t) is the 
fault-parallel shear stress draging on the base of the plate. For 
a homogeneous plate of shear modulus G, thickness averaging of the 
local stress-strain relation shows that 

d == G 3u/3y (2) 

A y 

where u(y,t) is the thickness average of the fault-parallel (and 
only nonvanishing) displacement. 

The shear drag is connected to the displacement and 
displacement rate, within the approximate treatment of coupling in 
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the model, through the linear Maxwellian relation 

(b/G)9i: x /9t + (h/Tp t x = 9u/9t - v o (y) (3) 

which describes the deformation of the viscoelastic asthenosphere . 
In (3), b is a short-time effective elastic coupling thickness. It 
is chosen as (7t/4) 2 H (Lehner et al, 1981) when G is the shear 
modulus of the elastic plate, in order for the final plate model to 
correctly respond to sudden stress release over the lithosphere 
thickness. Also, h/T] is a long-time viscous compliance, 
interpretable as a depth scale h over which asthenospheric shear 
takes place divided by asthenospheric viscosity T] . The velocity 
V Q (y) is the rate of motion of sub-asthenospheric material of the 
mantle. It is assumed to be a function of y, but independent of t, 
such that V Q (y) - V Q (-y) approaches the relative plate velocity V pl 
as y increases. The function V Q (y) represents the imposed steady 
driving motion although, as will be seen, the time dependence of 
deformation is remarkably independent of the detailed form of the 
function . 

On combining (1) , (2) and (3) , the model requires that 

«X + P9/9t) 9 2 u/9y 2 = 9u/9t - V o (y) (4) 

where oc = HGh/T) and P = bH ** (7tH/4) 2 . This is the equation governing 
the time and spatial distributions of deformation in the 
lithospheric plate, subject to a prescribed boundary condition at y 
= 0. (This boundary condition will be seen to be tied to V pl and 
the earthquake cycle repeat time T cy ) . 

At the plate boundary, y = 0 + , the thickness average stress G xy 
transmitted is assumed to be related to the net thickness average 
displacement 2 (u - D/2) through the stiffness k of the edge cracked 
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structure with freely slipping crack walls (Figure lb) . Here D = 
D (t ) represents the seismically accumulated slip of the presently 
locked shallow zone. The edge cracked structure in anti-plane 
strain may be regarded as a cross-sectional cut of the lithosphere 
at the plate boundary. For repeated earthquake cycles, D may be 
taken as a staircase function of time with periodic stepping. Fig. 
2. Thus, the appropriate boundary condition for (4) at y=0 + is 

G du/dy = 2k (u - D/2) . (5) 


The stiffness k is given by (see eq. (8) to follow) 


k = JCG/4H In [1/sin (JIL/2H) ] 


where L is the locked depth. This is consistent with the stiffness 
k given by Tse et al . (1985) for "line spring" modelling of 

partially locked plate margins. An appropriate choice of L would 
be the focal depth of earthquake ruptures. For the extreme case of 
through thickness rupture, L is equal to H and the stiffness k goes 
to infinity. Equation (5) then implies u = D/2 which recovers the 
simple boundary condition used, in an earlier analysis by Lehner 
and Li (1982) of strain accumulation in the earthquake cycle. 

To model periodic earthquakes with cycle time T cy and slip 
magnitude D^ the boundary slip D(t) may be decomposed into two 
parts, as shown in Figure 2. With the periodic saw-tooth function 
represented by a Fourier series, D(t) may be expressed as 

Oft 

D (t ) = (D 1 /7C) ^ (1/n) sin(2nJtt/T cy ) + [(1/2) + t/TJ (6) 

n=l 

In the steady state, it is assumed that the cumulative seismic slip 
agrees with the overall plate velocity V pl . Thus, D 1 = V pl T cy . The 
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second term of (6) may be understood as the value of 2u at y = 0 + 
for a particular solution for u in y>0 of form V pl t/2 + a function 
of y only, that reproduces the right hand side term V Q (y) of (4) . 
This solution has time rates representing steady rigid block motion 
and is otherwise of no interest in the present context. We are 
interested in the additional cyclically time dependent solution u 
to (4) with V Q = 0, in response to the first term of (6) . The 
solution of (4) subject to (5) may be accomplished by means of 
separation of variables, and the detailed procedure is given in 
Appendix A. The solution to within the function of y just 
mentioned is 


u(y,t) = (V pl /2) (t + T cv /2) 


cy' 


(7) 


CO 

+ V pl t r X <1/<W exP(-M„*> sin«O n t/t r - - NY) 


n-l 


where Y s y/Vp* 4y/JCH and the terms B n , \j/ n , M n and N n are functions 
of the circular frequencies G) n = 27Cnt r /T cy and of a parameter X 
(Appendix A) dependent on L/H. The Elsasser model relaxation time 
t r = P/a * (7t 2 H/16h) (T|/G) is essentially a fraction of the 
relaxation time for asthenosphere material; owing to the geometry 
of an elastic plate on a viscoelastic foundation, the time scale of 
lithosphere/asthenosphere coupling is appreciably longer (Rice, 
1980; Lehner et al, 1981) . All subsequent comparisons with data 
involve rates and hence are unaffected by the unwritten function of 
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The function u(y,t) as just calculated is the thickness 
averaged lithosphere displacement, and we must extract from it the 
displacement u s (y,t) at the Earth's surface. This is done as 
follows. Consider a plate ’ as in Fig. lb which has uniform shear 
modulus G and thickness H and has zero shear tractions on its base 
and on the vertical boundary below the locked ligament L. The 
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plate is subjected to a remotely applied shear force CTH, per unit 
distance along strike, and from solutions of the problem by Tada et 
al (1973), Turcotte and Spence (1974), and Tse et al (1985), this 
causes a uniform du/dy = c/G everywhere except at y = 0, where u is 
discontinuous due to slip below the locked ligament (i.e., on y = 
0, at depths > L) , such that 

u (y=0 + ) = -u (y=0“) = (2<JH/JtG) In [1/sin (7CL/2H) ] + D/2. (8) 

The last term reflects that the now locked ligament has previously 
slipped by amount D. Thus, the thickness averaged displacement 
increases linearly with y for y>0 and y<0, but has a jump 
discontinuity at the fault trace, y=0. The surface displacement u s 
mirrors this variation. Consider first the case D = 0. Then u s 
too varies linearly with y, and is equal to u, at large jy I, but u s 
is continuous at y=0. The jump discontinuity of u there is 
replaced by a continuous transition of u s over a length scale 
determined primarily by L, but by H also. This variation is 
described by the function 


S(y) = 1 


(Ity/2H) - In jsinh (;ty/2H) 


+ ^ 


sinh 2 (7Cy/2H) + sin 2 (7tL/2H) 


(9) 


In [1/sin (JIL/2H) ] 


which can be extracted from the references above, and varies from 
S=0 at y=0 to S=1 at large y, such that (for y>0) 

u s (y) = .[u(y) - u(0 + )] + u (0 + ) S(y) (10) 

= u(y) - u(0 + ) [1 - S (y) ] 

When there is previous slip by amount D the equation remains 
correct if each u s and u, for y>0, has subtracted from it a rigid 
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motion D/2 so that there results 

u s (y) = D/2 + [u (y) - u(0 + )] + [u(0 + ) - D/2] S (y) . (11) 

As examples, when L = H/4 one has S (O) = O, S (L) = 0.54 and S (H) = 
0.96, and when L = H/2, S(O) - O, S(L) = 0.72 and S (H) = 0.94, 
illustrating the rapid approach to unity. 

For the case when the lithospheric plate is coupled to a 
viscoelastic asthenosphere, we can use the foregoing analysis. as a 
basis for estimating u s (y,t). That is, we assume as in the exact 
analysis reflected by eqs . (10, 11) above that u s is essentially 

equal to u, except that the jump discontinuity in u at the fault 
trace shows up as a gradual variation in u s spread out by the same 
function S (y) . Thus we calculate u g from 

u s (y,t) = D(t)/2 + [u (y, t) - u(0 + ,t)] 

+ [u(0\t) - D(t)/2] S (y) (12) 

where D(t) is given as in Figure 2 and u(y,t), and thus u(0 + ,t), by 
the solution in eq. (7) . Rates of fault parallel surface 
displacement and shear strain are then 

u s = 3u s (y,t)/3t , Y s = 3 2 u s (y,t) /3ydt (13) 

and these quantities are shown in various subsequent plots. 

Since u(0 + ) - D/2 is proportional, by eqs. (5) and (8), to the 
net force 0(0) H transmitted across the locked ligament, eq. (12) 
expresses our assumption that it is this net force, together with 
geometric dimensions L and H as included in S(y), which determines 
the perturbation of the surface displacement and strain profile due 
to the freely slipping deep fault surface below. 
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For the present, we observe that the solution derived here has 
the form 

u s (y,t) = (V pl /2) F 1 (y/H, t/T , L/H, t r /T ) 

(14) 

Y,(y/t) = <V pl /H) F 2 (y/H, t/T cy , L/H, t r /T cy ) 

where F x and F 2 are universal dimensionless functions of their 
dimensionless arguments. Further as y — » °o, F x -» 1 and F 2 -> 0. 
Also, Fj^ is asymmetric and F 2 symmetric in y. These forms will 
enable interested readers to translate specific results given in 
later plots to some other sets of parameters. 

The governing equation (4) of the deformation of the plate is 
based on the modified Elsasser model (Rice, 1980; Lehner et al, 
1981) . A characteristic of this plane stress model is in the use 
of thickness averaged quantities. Based on through thickness 
ruptures, Lehner and Li (1982) computed displacement profiles using 
the modified Elsasser approximations. Their results were compared 
with finite element calculations by Cohen and Kramer (1984), who 
suggested that the thickness averaging of parameters contributed a 
major source of inaccuracy in the prediction of the surface 
displacement field. 

In the present model, the computed surface displacement (and 
subsequently the displacement rates and strain rates) takes into 
account the shallow locked depth L<H, which is embodied in the S (y) 
terms in (10) - (12) . We expect that the surface deformation field 
would then be most sensitive to the effect of continual shear below 
the locked ligament so that the above objection to the approach 
does not apply. We demonstrate that such is the case by 
considering (Appendix B) the extreme model of plates loaded at 
their base by uniform imposed displacements of opposite sign for 
y>0 and y<0, consistent with a uniform plate velocity V 1 . (This 
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is a severe test since the base traction rates X x are then 
localized near the fault trace, whereas our solutions show much 
broader loading zones (hence conditions more consistent with our 
use of the solution for the remotely loaded edge-cracked strip) 
well into the interseismic period; see the final figure in Appendix 
B.) 


This boundary value problem has been solved exactly by 
conformal mapping techniques (Li and Lim, unpublished manuscript, 
1986) and the relevant expressions for displacement rate and stress 
rate are given in Appendix B. For this model, eqn. (B3) and eqn. 
(B5) give the exact form of the surface displacement rate profile 
u g (y) and the surface strain rate profile Y s (y) . If we now 
calculate for the base loading the associated thickness average 
quantities, eq. (12) and (13) then provide approximate estimates of 
the surface deformation consistent with our present procedure. 

Comparisons of the exact u s (y) and Y s (y) with the approximate 
solutions, shown in Appendix B, reveal that the surface deformation 
fields calculated with (12) and (13) would be reasonably accurate 
in the case where the locked zone is shallow. We have made further 
comparisons between the exact and approximate models for the shear 
strain rates at the fault trace, and found that Y exact /Y ap prox = 0.82 
- 0.91 for the locked zone depth and the lithosphere thickness of 
the best fit model for the San Andreas to be discussed shortly. 
This exercise demonstrates that the surface deformation field near 
the fault trace is more sensitive to the locked depth and the net 
driving load at the boundary than to the details of the coupling to 
the asthenospheric foundation. Thus, the approximations made in 
using thickness average displacement to estimate surface 
deformation by means of (12) and (13) is not expected to produce 
intolerable inaccuracy. (We assume that the model provides a 
reasonable account of time dependent force transmission at the 
margin and of thickness averaged displacement beyond the short term 
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post seismic period). Note that Cohen & Kramer's comparison between 
their FE calculation and the modified Elsasser formulation is based 
on the full thickness rupture L = H analysis due to Lehner and Li, 
and is shown to be the worst case for the Yexact/Yapprox ratio - 
Further, the loading mode discussed provides a severe test of our 
procedures, which we know to become exact when base loads are 
instead applied remotely from the cracked section. 

An unsatisfying aspect of the base loaded model shown in Fig. 
B1 is in the uniform imposed displacement (or velocity) , implying a 
rigid foundation. This unrealistic feature, however, is overcome 
in the present Elsasser formulation. 

It may be of interest to note that the 
lithosphere/asthenosphere coupling model of Savage and Prescott 
(1978) and its derivatives (Cohen and Kramer, 1984; Thatcher, 1983) 
is based on elastic plates overlying a viscoelastic half space in 
which the deformation in the shear zone below the seismogenic depth 
within the elastic plate is imposed kinematically as spatially 
uniform slip at a constant rate consistent with the remote plate 
velocity V pl . This imposed uniform shear slip distribution, rather 
than one determined as here by coupling between the lithosphere and 
asthenosphere, strongly moderates the time-dependence of surface 
deformation. An implication is that for shallow faults, L«H, 
coupling is predicted to be almost non-existent in their model 
since the near fault surface strain field would be almost 
completely controlled by the imposed dislocation rate. Indeed, 
Savage and Prescott (1978) reported that in their model, "the 
effect of asthenosphere relaxation is important only if the depth 
of the seismic zone is comparable to the thickness of the 
lithosphere". In contrast, the present model predicts response of 
the deep aseismic shear zone and it is found to respond in a 
spatially (in z) and time varying fashion as shown in Fig. 4 (see 
further explanation below) . 
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Theoret i ca l — re.sult S , — model parame ters, and comparisons with 

time-dependent strain data from California 

The complete surface strain rate at various fractions of an 
earthquake cycle is shown as a function of distance from the fault 
in Figure 3a, b. These figures compare results for two extremely 
different models, (a) locked through the entire depth H and (b) 
locked only over 0.3H. Immediately after the last rupture, the 
strain rate is high and localized near the fault due to rapid 
relaxation of the asthenosphere . It decays in amplitude and 
spreads spatially with increasing time. The rate of spreading is 
associated with the relaxation time of the asthenosphere. 
Observational support of such spreading of deformation is suggested 
by Thatcher (1983) on the basis of comparing the northern and 
southern sections of the San Andreas fault, and his composite 
estimate of strain rates at the fault trace for the two regions 
suggests decay in time. As expected, comparison between Figures 3a 
and 3b shows that a shallow locked segment tends to cause more 
localized deformation (at the fault) than one with most of the 
fault depth locked during the earthquake cycle . 

To better understand the localized time-dependent straining at 
the fault, we recall that the relaxation of the asthenosphere 
causes a reloading of the plate boundary. This reloading causes a 
gradual shear flow of the fault zone at depths below the locked 
ligament. To illustrate, we show in Figure 4 the crack face slip 
(measured arbitrarily from zero immediately after an earthquake) at 
various fractions of a cycle time for the case in Fig. 3b. Note 
that the rate of slip as well as the surface strain rate slows down 
as the asthenosphere relaxes . For the case illustrated more than 
70% of the total postseismic slip has occurred in the first half of 
the cycle. This is consistent with the post-seismic high strain 
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rate shortly after the earthquake rupture as mentioned earlier. 

We wish to make our description of the earthquake loading 
process, which (within the framework of the model adopted) means 
our choice of the five parameters V pl , L, H, T cy and t consistent 
with seismic, geologic and geodetic evidence from the San Andreas 
fault region. Our approach is to choose V pl , L and T cy on the basis 
of geologic and seismic considerations and then to choose the less 
well constrained parameters H and t r to fit the data set assembled 
by Thatcher (1983) for time dependent post-earthquake straining at 
a variety of locations along the San Andreas fault. Predictions 
based on the resulting parameter set are then compared to 
comtemporary fault-parallel displacement rates as a function of 
distance from the fault trace in a few places. 

In the case of V pi' it is accepted that the overall velocity 
between the Pacific and North American plates is approximately 55 
mm/yr (Minster and Jordan, 1978) . However, evidence from long term 
fault offsets summarized recently by Minster and Jordan (1984) and 
Stuart (1984/85), based on work by Sieh and Jahns (1984) and Weldon 
(1984) , suggests that only 35 mm/yr should be regarded as taken up 
by the San Andreas fault, and we therefore set V pl = 35 mm/yr. 

In northern California great earthquakes like the 1906 rupture 
have been estimated to have a repeat time from 70 to more than 180 
yr (Savage, 1983) . A longer repeat time of 225 yr. based on 
trenching data was obtained by Hall et al. (1982) . In southern 
California, data from extensive trenching at Pallet Creek (Sieh, 
1984) leads to an average recurrence interval between 145 and 200 
yr for large earthquakes, based on excavated evidence of 12 
earthquakes between approximately 260 and 1857 A . D . Evidence from 
other locations along the 1857 rupture zone (Sieh and Jahns, 1984; 
Stuart 1984/85) suggests the possibility of a variable average 
recurrence interval from one segment to another, ranging from 
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approximately 100 to 300 years. We have fixed T cy = 160 yr, since 
this value is reasonable in terms of observations and since 
variations by ± 25% (expanding the range to 120 to 200 yr) have 
only very modest effects on the results which we show. 

In terms of the model of a locked shallow portion of the crust 
with accumulating aseismic slip below, it seems logically 
consistent that the nucleation depth for a large earthquake should 
be identified with the highly stressed region at the border of the 
slipping and locked zones, i.e., with the crack tip location. This 
would argue for choosing L in the range of 8 to 10 km on the basis 
of seismically determined nucleation depths. For example, the 1979 
Imperial Valley earthquake is inferred to have initiated at about 8 
km depth (Archuletta, 1982), the 1966 Parkfield earthquake at about 

9 km depth (Lindh and Boore, 1981), and both the 1979 Coyote Lake 
and 1984 Morgan Hill earthquakes along the Calaveras fault at about 

10 km depth (Lee et al., 1979; Bouchon, 1982; Bakun et al., 1984). 
Also, Thatcher (1975) inferred a focal depth of not greater than 
about 10 km for the 1906 San Francisco earthquake based on analysis 
of geodetic data. 

The nucleation depths discussed are towards the bottom of the 
seismogenic zone along the San Andreas, i.e., the zone extending 
typically to 12 to 14 km depth from which seismic signals emanate 
(Eaton et al., 1970; Wesson et al., 1973; Sibson, 1982). 

In order to associate L with representative nucleation depths, 
we take L = 9 km as one case for more detailed study. However, in 
our present simplified modelling the locked zone depth L is 
regarded as fixed throughout the earthquake cycle. In reality, L 
will be time dependent if the slipping region gradually penetrates 
into the locked zone throughout the cycle. In that case the 
earthquake nucleation depth will underestimate a representative 
value of L for, say, mid-cycle straining conditions. We suggest 
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that a choice of L = 11 km, which we take as a second case for 
detailed study, reasonably represents a locked zone depth in the 
presence of upward penetration. 

As a first perspective on this choice, we note that the shear 
crack-like configuration implies that the locked region ahead of 
the crack tip is of mechanical properties such that it tends to 
resist deformation, but nevertheless ultimately breaks, in that it 
slips subsequently at reduced strength, when appropriate local 
stress levels are reached. Such properties define a deformation 
weakening material, i.e., a potentially brittle material, and are 
plausibly associated with material of the seismogenic zone. Thus 
the greatest plausible depth of locked zone would correspond to the 
12 to 14 km base of the seismogenic zone; such would give the 
locked zone depth early in the cycle, whereas late in the cycle the 
depth has diminished to the 8 to 10 km nucleation depth range. 
These considerations suggest an upward motion of the locked zone on 
the order of 4 km throughout the earthquake cycle, and thus a mean 
depth around 10 to 12 km (we take L = 11 km) , if indeed the locked 
zone depth does start off the cycle approximately coincident with 
the maximum seismogenic depth. 

An independent but remarkably supportive perspective is 
provided by theoretical instability modelling. The shear crack and 
slip weakening instability models mentioned earlier all show a 
gradual upward penetration of the slipping zone. The recent work 
by Tse and Rice (1986) is perhaps the most realistic of that series 
in adopting a slip rate and slip history dependent form for the 
fault surface constitutive law based on laboratory experiments, and 
in using what (very incomplete) data is available on temperature 
and normal stress dependence of constitutive parameters, together 
with a San Andreas fault geotherm, to provide a depth variable 
description of material response. This (forward) modelling implied 
a transition from inherently stable to potentially unstable slip at 
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around 11 km depth. Further, the results for slip histories 
throughout whole earthquake cycles calculated by Tse and Rice 
(1986) show clearly the development of a shallow locked zone which 
is gradually penetrated by slip below. Their simulations allowed 
for a broad range of critical slip weakening distances and 
different options for describing uncertainties of constitutive 
response; all were done with V pl = 35 mm/yr and the resulting 
values for T cy were between 50 and 160 yr. Examining simulations 
showing T cy > 80 yr, one observes that the slipping zone below 
reached up to around 8 to 9 km depth when the deep slip reached 25% 
of V pi T cy 25% of its full cycle value), moved up to 7 to 8 km 
at 50% of V pl T cy , then to 6.5 to 7.5 km depth at 75% of V pl T cy , and 
to around 6 km just before instability. We may summarize by saying 
that the Tse— Rice modelling had an 11 km deep seismogenic zone and 
a locked depth at nucleation around 6 km (both a few km shallower 
than estimated above) and a locked zone depth around mid— cycle on 
the order 1 to 2 km deeper than that at instability, and 3 to 4 km 
shallower than the maximum seismogenic depth. These modelling 
results are not inconsistent with the simple description given 
above, and suggest that a representative locked zone depth at L = 
11 km, 2 km below a typical nucleation depth, may account 
approximately for upward penetration. A feature not accounted for 
is that upward penetration would cause the strain near the fault 
trace to slightly increase even if the net force transmitted across 
the plate boundary were not increasing with time. This effect may 
act to slightly counteract the general trend of decreasing strain 
rate throughout the cycle that our present modelling implies . 

For definiteness, we choose two cases L = 9 km and L = 11 km 
in what follows along with V pl = 35 mm/yr and T cy = 160 yr. There 
are no comparably good apriori constraints on the appropriate 
values of lithosphere thickness H and viscoelastic relaxation time 
t r for purposes of describing the earthquake stressing process . 
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Thus as a general constraint on these model parameters for the 
San Andreas fault as a whole, Figure 5, we first compare the 
surface strain rate Y s computed at the fault with Thatcher's (1983) 
compilation of data from Pt. Arena, Fort Ross, Pt . Reyes, Shelter 
Cove, S.F. Peninsula, Carrizo Plain, San Fernando-Mo jave, Palmdale 
and Coachella Valley. Locations are shown in Figure 5, and 
Thatcher's plot of the data is shown in Figure 6 with his estimated 
error bars indicated. In reference to figure 6a, taking the locked 
zone as 9km, we found by a series of runs with different parameter 
values that a good fit to the data is given by t r =12 yr and 
H=22.5km. The results of our calculation in that case are shown by 
the thick solid line, which we term our "reference" curve. To 
explain why fitting the data serves to constrain both t r and H, we 
now examine the alteration to this reference case by varying one at 
a time. For the case of lowering the relaxation time to 6 yr, it is 
seen that the curve rotates by increasing the strain rate 
significantly at short time after the last rupture and drops below 
the reference curve approximately 20 years after the last rupture. 
While the short time prediction can arguably fit the geodetic data, 
the long time curve appears to underestimate the strain rate, lying 
below every data point except PR. As we shall see, the assumption 
of a thinner lithospheric thickness shifts the curve up at all 
times, and thus may help to improve the t r =6 yr curve fit at long 
times, but would then aggravate the short time prediction. A 
similar analysis done for a much higher relaxation time of 18 yr. 
has the effect of rotating the curve downwards at short times and 
upwards at long time. While at short times, it passes through the 
error bars of the PA and FR points which may suggest a suitable 
fit, the long time prediction appears to over-estimate the strain 
rate data, with a clear exception of CP. One may again attempt to 
salvage this fit by using a larger H, which would indeed help in 
the long time, but would then aggravate the short time fit. The 
above attempts to study the effects of varying t r suggests that t r 
should be bounded by a lower limit of 6 yr. and an upper limit of 
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18 yr. The best fits appear to be in the range of 10 to 14 yr. 

A study of the effect of varying the plate thickness from the 
reference value of 22.5km was also carried out. These results are 
also plotted in figure 6a and show an upper value of 30km (which 
puts the prediction on the low side compared to data) and a lower 
value of 15km (which puts the prediction on the high side) . Again 
one may attempt to improve these fits by changing the relaxation 
time, but this always results in improving one part of the curve at 
the cost of making the other part much worse because of the 
rotating effect as discussed earlier. Thus the fitting of the data 
seems to constrain both t r and H; significant alternations of one 
of the two from the value for our good fit reference curve cannot 
be compensated for by alteration of the other. 

We have also analyzed the case for the locked zone L=llkm. 
The results show a good fit with t r =14km and H=20km, shown as the 
thick solid reference curve in figure 6b. In this figure we have 
also shown two fits of varying H up to 27.5km and down to 12.5, 
which somewhat bracket the data points . Again, we conclude that 
large variations of t r or H from the reference values cannot give 
good fits for the reasons given in the previous paragraphs. 

Since displacement rate profile predictions will be performed 
for various selected locations in this paper, it would be of 
interest to gain some insight into the influence of the model 
parameters on the predicted profiles. For this purpose, the same 
set of parameters used for the strain rate decay study in figure 6a 
was employed to compute the corresponding surface displacement rate 
profiles for a typical 'contemporary' time, t = 100 yr, and this 
result is shown in figure 7a. It shows that the displacement rate 
generally increases (in absolute values) with a higher relaxation 
time t r and a smaller lithospheric thickness H, at least for this 
late (more than 60% of) cycle time. The variation of locked zone 
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depth may also affect the displacement rate profiles, although the 
effect is ralatively small by varying L=6km to L=12km, as is shown 
in figure 7b. 

Combination of eqs . (1) and (2) shows that the curvature 
3 2 u/3y 2 of a plot of u versus y is proportional to the rate of 
change T x of mantle drag stress on the base of the plate. That is, 
it is proportional to the local rate of loading of the plate. The 
plots in Fig. 7 show u g rather than u but the two are virtually 
indistinguishable for y > H. The fact that the plots in Fig. 7 
still show appreciable curvature out to distances of order 3H (the 
plots all asymptote to a common level V pl /2 = 17.5 mm/yr as y — » °°) 
show that the region of base loading is very broad. This is in 
sharp distiction to the case discussed in Appendix B and, we argue, 
justifies our approximations embodied in eqs. (5) and (12) in 
analysis of the edge cracked strip. 

Let us summarize as follows. If on the basis of typical 
rupture nucleation depths we assume the locked zone depth L = 9 km, 
and take V pl = 35 mm/yr and T cy = 160 yr, then the data of Fig. 6 on 
time dependent straining throughout the earthquake cycle is 
sensibly fit by assuming a lithosphere thickness H = 20 to 25 km 
and a relaxation time t r = 10 to 14 yr. Substantial deviations of 
H and t r from these ranges were shown to not allow acceptable fits 
to the data. Similarly, assuming L = 11 km to account 
approximately for the expected upward motion of the base of the 
locked zone during the cycle, i.e., an effective L greater than a 
typical nucleation depth, the data is sensibly fit by assuming H = 
17 to 23 km and t r = 12 to 16 yr. 
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Discussion : li thosphere thick ness, possible crustal asthenosohere 
and relaxation time 

It is useful to compare the inferred H and t r values to some 
constraints obtained from other observations. This provides some 
confidence in our model predictions and sheds light on the larger 
tectonic implications of our results . 

A wide range of lithosphere thickness values have been assumed 
in previous studies. Averages of oceanic and continental 
lithosphere thicknesses suggested by Stacey (1977) come to about 90 
km, whereas Anderson (1975) used a 50 km oceanic lithosphere and 
Lehner et al. (1981) used a 75 km average lithosphere for generic 
illustration of viscoelastic coupling effects within the Elsasser 
framework. In a recent re-evaluation of isostatic rebound at Lake 
Bonneville, Nakiboglu and Lambeck (1983) used 28 to 30 km. 

At the other extreme, and specifically for a lithospheric 
thickness along the San Andreas fault, Thatcher (1983) suggests 
that the elastic lithosphere may have a thickness of only 10 km. 
This was motivated by an attempt to rationalize the same data as 
in Figure 6, although our work here based on what we think to be a 
physically more valid model, without artifical kinematic imposition 
of deep fault slip, shows clearly that the data is compatible with 
much thicker values, 17 to 25 km. Further, Turcotte et al. (1984) 
pointed out that Thatcher's thin lithosphere would require a 
significantly higher heat flow than the relatively high value found 
over a broad region near the San Andreas fault, and they suggested 
a 30 km thick three layer lithosphere model which they coupled to 
the deeper surroundings by an Elsasser approximation. Their middle 
layer was intended to represent an intracrustal asthenosphere . 
However, a deficiency of their modelling is that they included no 
account of ongoing aseismic slip on the downward continuation of 
the locked seismic zone. Thus, while Turcotte et al. (1984) showed 
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that they could choose parameters of their layer model to fit the 
sort of time-dependence of strain near the fault trace shown in 
Figure 6, we expect that in the absence of the strain concentrating 
effect of deep slip at the plate margin their model would not 
describe adequately the variation of surface strain or displacement 
rate with distance from the fault trace. They did not show results 
for what their model implied about such variations . 

Finally, Tse et al (1985) have modelled geodetic data from the 
Parkfield region and adjacent central California creeping segment 
of the San Andreas fault by a model similar to what we examine 
here, involving elastic plates locked over part of the depth range 
at the plate margin, but freely slipping elsewhere. They assume a 
stress free base for the plates and thus load the system by 
remotely applied forces, as in the Turcotte and Spence (1974) crack 
modelling, but allow the geometry of the locked ligament to vary 
along strike and to not exist at all within the central region of 
the creeping zone. They show that the fit to data is better with H 
in the 30 to 40 km range than for larger values, but did not 
examine smaller values. Thus we can interpret their work as 
suggesting that H is less than approximately 40 km. 

Evidently, the lithosphere thickness values of order 17 to 25 
km that we infer are not incompatible with other studies. Further, 
these thicknesses are comparable to or less than representative 
crustal thicknesses, in that Oppenheimer and Eaton (1984) infer a 
Moho depth slightly above 24 km for the northern San Andreas 
between the PA (Point Arena) and CP (Carrizo Plain) symbols in 
Figure 5, with a thickening towards 26 km within approximately 50 
km of those two points and a modest thinning towards 23 km at PR 
(Point Reyes) . The tectonic significance of this observation is 
that the viscoelastic asthenosphere which enters our modelling must 
be presumed to correspond to the lower portions of the crust and 
perhaps a crust-mantle transition zone. 
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There may, of course, also exist beneath the San Andreas fault 
region a classical deep asthenosphere lying beneath a more elastic 
layer of upper mantle material, but we tentatively conclude that it 
is a lower crustal asthenosphere rather than that intra-mantle 
asthenosphere which contributes to time variations in the crustal 
deformation rate throughout the San Andreas earthquake cycle. We 
remark now and will return to the observation from Oppenheimer and 
Eaton (1984) that the crustal thickness decreases towards the ocean 
(i.e., southwest) side of the San Andreas fault, and increases 
toward the continental side. Near Point Reyes, this occurs at a 
steep rate of about 1.2 km change in crustal thickness per 10 km 
distance perpendicular to the fault (but they show thickness 
contour lines only on the continental side there and indicate that 
the uncertainties are greater than elsewhere) , and more typically 
at a rate of 0.6 to 0.8 km change in crustal thickness per 10 km 
perpendicular distance. One infers, spurred on by some puzzling 
geodetic observations of Prescott and Yu (1986) to be discussed 
shortly, that to the ocean side of the San Andreas near Point 
Reyes, the already thin crust may possibly become too thin (hence 
too cool at its base) to sustain a crustal asthenospheric layer 
near its bottom. 

An initially unexpected consequence of the inferred 
asthenosphere location in the lower crust is that the Elsasser type 
approximation for the viscoelastic coupling becomes yet more 
suitable since the relaxation would then seem to be confined to a 
relatively thin channel. 

The relaxation time t r is proportional to HT|/hG, and there is 
considerable uncertainty associated with an apriori choice of a 
suitable asthenospheric viscosity T] and depth scale h over which 
velocity gradients occur in viscous relaxation. Previous attempts 
to estimate T| have generally assumed that the asthenosphere 
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involved is of the classical intra-mantle type and are of uncertain 
relevance if a thin lower crustal asthenosphere is the correct 
picture. Nevertheless, with regard to T|, Nur and Mavko (1974) 
inferred 5 x 10 18 Pa-s from post seismic relaxation whereas 
Thatcher and Rundle (1979) and Thatcher et al (1980) suggested 
values of 4 x 10 19 Pa-s and 1 x 10 19 Pa-s, respectively, for 
post-seismic lithosphere rebound at underthrust zones. Walcott 
(1973) inferred 5 x 10 19 Pa-s from isostatic rebound data at Lake 
Bonneville, and a recent reevaluation of that case by Nakiboglu and 
Lambeck (1983) suggested T] = 1.5 to 3.4 x 10 19 Pa-s when a 28 to 30 
km elastic plate is assumed to lie over a uniform viscoelastic 
half-space, and 2.1 to 5.8 x 10 18 Pa-s when the plate lies over a 
100 km thick channel asthenosphere. 

The crustal shear modulus near the San Andreas is G = 35 GPa, 
based on the average p-wave velocity of 6.0 km/sec (Oppenheimer and 
Eaton, 1984) and assuming a Poisson ratio of 0.25 and specific 
gravity of 2.9. Thus, the relaxation time t r = 13 yr (middle of 10 
to 16 yr range) implies T| = 2.3 (h/H) x 10 19 Pa-s. We have 
inferred that H = 17 to 25 km. Using H = 21 km and, since the 
asthenosphere is assumed to be confined to the lower crust and 
crust-mantle transition zone, assuming that h may range from 2 to 
10 km, we conclude that viscosities of 2 x 10 18 Pa-s to 10 19 Pa-s 
are implied. These are in a range compatible with other estimates 
as just summarized. We have no independent check on their 
suitability or, indeed, on the adequacy of a linear 
characterization of viscoelastic response. 


Contemporary fault parallel strain and displacement rates 

Here we compare our model predictions, based on the parameters 
already chosen, to data other than that used to constrain the 
parameters. The data consists of geodetically measured 
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contemporary fault parallel surface displacement rate profiles as a 
function of distance from the San Andreas fault, at three locations 
along the fault, and contemporary strain rates averaged over 
variable-sized areas at one location. Within the context of a 
model that involves uniform conditions along strike, and symmetry 
relative to the single fault strand representing the San Andreas in 
our surface plate, we believe that the comparison of prediction to 
data here is supportive of our model. Nevertheless, the 
displacement rate profiles at particular locations along strike, 
constructed from data for individual geodetic markers, are strongly 
affected by local conditions. Sometimes these suggest pronounced 
asymmetry of material properties or geometry relative to the San 
Andreas and show the effects of adjacent fault strands. Also, the 
data include features that suggest large non— uniformities along 
strike or, perhaps instead, inaccuracies relating either to marker 
instability relative to the crust below or to systematic 
measurement error. In one case, we note changes in model 
parameters which would improve the fit, but this is not always 
feasible within the symmetry of our model relative to the San 
Andreas fault . 

Palmdale area (point PD in figure 5 and figure 8 and 9) : 

The displacement rate (average velocities, 1973-1983) data 
was deduced by King and Savage (1984) for 11 of the stations in 
Figure 8, from line length measurements using the "outer 
coordinate" solution method which minimizes the rms displacement 
normal to the San Andreas fault. Based on the line length changes, 
King and Savage concluded that there was no clear evidence of 
surface fault slip in the Tehachapi region. This data set has a 
relatively large scatter, and is plotted in figure 8 with the 
displacement rate zeroed on the San Andreas fault. The two 
theoretical curves shown are based on the parameters of the best 
fits to the strain rate decay data set shown in figure 6a and 6b. 
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They are computed at t = 123 yr (1857 + 123 = 1980), within the 
time period when the time measurements were made. Within the 
scatter of data, these parameters give apparently reasonable 
predictions to displacement rate profiles in the Palmdale area. 
The data is probably affected by the presence of the Garlock fault . 
King and Savage fit it with a buried screw dislocation model and 
conclude that a slightly better fit is obtained if they assume deep 
slip on the Garlock fault (at slightly less than half the rate they 
infer on the San Andreas), rather than slip on the San Andreas 
only. Their inferred locking depth is a poorly resolved 15 to 20 
km, which we regard as questionable on seismic grounds and, in 
light of our results with much shallower locking depths, probably 
due to shortcomings of the buried dislocation model which we have 
discussed. 

The strain rate averaged over various sub-nets of the network 
shown in Figure 8 (see figures 5 and 6 of King and Savage for 
locations of the sub-nets) was calculated by King and Savage (1984) 
and is listed in Table 1. Also listed is data for the Los Padres 
geodetic network, lying just to the NW of the region shown in Fig. 
8 and mostly to the ocean side of the San Andreas (Savage, 1983), 
as summarized by McGarr et al. (1982) . We have calculated the 
corresponding strain rate from our model using the same two sets of 
parameters mentioned earlier. The calculation was done based on 
the centroidal distance y of the net from the San Andreas fault 
(positive on continental side) and the root mean square distance. 
(Ay) rms of the stations in the net, as given by McGarr et al . 
(1982), and these values are also given in Table 1. The strain 
rate calculations were based on a procedure attributed to Savage by 
McGarr et al., using 

Y = [u s (y + (Ay) rms ) - u s (y - (Ay) rms )]/2 (Ay) 


(15) 
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They implemented this formula with the buried screw dislocation 
model . 

Our results based on our sets of reference parameters are 
tabulated as y(a) for the L=9km case and y(b) for the L=llkm case. 
These calculated results are seen to be on the high side when 
compared to the observed strain rate. However, we have not 
attempted to fit this set of strain data, and the comparison can 
be improved by varying the model parameters . Owing to the 
uncertain effect of the Garlock fault, it is not clear that a 
fitting is justified. Nevertheless, Figures 7a and b show how 
slopes of the displacement rate profiles, i.e., the strain rates, 
can be altered by variation of L, H, and t r , and, of course, the 
curves shown all have amplitude directly proportional to V pl . We 
show as y(c) in Table 1 the result which correspond to V pl = 32 
mm/yr, L = 9 km, H = 25 km, t r = 12 yr, and T cy = 160 yr just to 
emphasize that a good fit is possible with a locked depth that, we 
have argued, can be justified on seismological and material 
property considerations. In comparison, McGarr et al (1982) find 
that with the buried screw dislocation model the data implies a 
locked depth, again poorly constrained, around 22 km. 

We may note also that the buried screw dislocation fits of 
King and Savage and McGarr et al . imply contemporary deep slip at 
rates of order 20 to 25 mm/yr. Our modelling shows that while 
points at the base of the surface plates, far removed from the 
fault trace, have relative motion at rates of V pl = 35 mm/yr, the 
average contemporary slip rate on the freely slipping deeper 
portion of our fault zone, 2 [H/ (H-L) ] u (0 + , t) , is only of order 4.9 
to 6.4 mm/yr. The reason for this great difference will be evident 
from figures 3b, 4 and 7: at large times since the last earthquake 
the average deep slip rate diminishes considerably, but the 
accommodating region of the asthenosphere broadens out to several 
times the plate thickness. 
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Point Reyes area (point PR in Figure 5 and Figures 10 and 11) : 

Data for the period 1972-1982 from geodetic networks at Point 
Reyes (38.1°N), Santa Rosa (38.3°N) and Napa (38.0°N) was analyzed 
by Prescott and Yu (1986), updating data given previously by 
Prescott et al (1979). These networks are shown in Figure 10. 
Because of the proximity in latitude of the three mentioned 
locations, the displacement rate values from these networks provide 
a profile of fault parallel rate as a function of distance from the 
San Andreas. Figure 11 shows the component of surface displacement 
rate parallel to the fault trace, as compared to data assembled for 
all stations of each network by Prescott and Yu (1986) . Their 
displacement data is analyzed as in Prescott (1981) , in that a 
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rigid body rotation has been chosen so as to mimimize displacements 
perpendicular to the fault trace. We have added a rigid 
translation to the data of Prescott and Yu (1986) so that zero 
fault-parallel displacement occurs at the San Andreas. The 
triangles are from the Geyser network to the north (Figure 10) . 

Our theoretical curves shown in Figure 11 with this data use 
the same sets of best fit parameters as deduced from the strain 
rate decay fit (figure 6a, b heavy solid curves) . That is, the 
parameters are chosen independently of the Pt . Reyes area data set. 
The profiles are computed at t=73 yr. (1906+73=1979), within the 
time period when the line length measurements were made. The 
predicted profiles appear to match the observed displacement rates 
very well on the NE side of the fault up to distances beyond the 
Rogers Creek fault extending to y = 40 km. Beyond the West Napa 
fault, however, the data show clear deviation with higher rates 
than the predicted profiles. This may reflect accommodation of the 
part of the 55 mm/yr overall plate motion not attributable to the 
San Andreas. Similarly, while the predicted profiles match the two 
data points close and to the SW of the San Andreas, the other data 
points corresponding to the Point Reyes Head and Farallon Islands 
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station velocities show marked deviation from an asymmetric 
displacement rate field about the San Andreas. We cannot judge the 
reliability of the data. As Prescott and Yu (1986) note, if valid, 
they suggest very much smaller shear strain rates to the SW of the 
fault than to the NE . 

We have observed earlier that our relaxing asthenosphere is 
inferred, by fits to the strain rate decay data, to lie in the 
lower crust. Further, the seismically inferred crustal thickness 
by Oppenheimer and Eaton (1984) suggest that the crust is already 
slightly thinner near Point Reyes than elsewhere along the northern 
San Andreas fault and seems to be diminishing in thickness at an 
unusually high rate as one moves towards the fault from the NE 
side. Thus, the possibility arises that the crust becomes 
sufficiently thin on the SW side near Point Reyes that its (now 
shallow) lower regions are too cold to allow viscous relaxation. 
This line of argument suggests that the crustal asthenosphere may 
get pinched off to the SW side of the fault, so that the effective 
lithospheric thickness becomes much greater on that side and 
whatever deeper relaxation occurs there takes place in a classical 
sort of intra-mantle asthenosphere. This should result in a 
pronounced asymmetry of surface straining relative to the fault 
trace, but we do not yet know how to model it to compare against 
the data. Also, the discussion is speculative and driven by the 
two data points (Farallon Islands, Point Reyes Head) of uncertain 
significance. For example, aseismic undersea fault regions cannot 
be ruled out and could influence either data point. 

Prescott and Yu (1986) develop an elaborate kinematical model 
to fit the data based on a distribution of buried screw 
dislocations in a half space. Their preferred distribution (model 
E in their paper) involves a pair of oppositely slipping screw 
dislocations at 6 and 10 km below the San Andreas, so as to 
accomplish spatially uniform slip at 10 mm/yr over the 6 to 10 km 
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depth range, and to this they add an infinite number of 
inf initessimally slipping screw dislocations (i.e., a continuous 
distribution) at 10 km depth, spread out uniformly over the 50 km 
lateral distance from .the San Andreas to the West Napa faults, so 
as to sum to 30 mm/yr total motion. In comparison, our theoretical 
predictions as shown in Fig. 11 suggest that the observed broad 
deformation profile between the San Andreas and West Napa fault may 
reasonably be interpreted as the result of ongoing mantle motion, 
together with laterally spreading viscoelastic asthenospheric 
relaxation from the great 1906 earthquake, taking place below an 
elastic plate of order 20 km depth. These features appear 
naturally in our physical (as opposed to kinematical) modelling, 
based on parameters chosen from other data sets. 

Coachella Valley - Salton Sea Area (point CV in Figure 5 and 
Figures 12 and 13) : 

These networks cover a rather complicated area, with the San 
Jacinto fault and the Elsinore fault running parallel to the San 
Andreas fault, see Figure 12. These lines have been analyzed by 
Savage (1983) , who obtained displacement rates (averaged over 
1972-1981) as shown in Figure 13. The computed profiles are based 
on the same set of parameters which best fit the strain-rate decay 
data for the whole San Andreas fault (figure 6a, b) . The profiles 
are calculated at 94% of cycle time, assuming the region has a 
relaxed foundation. However, the earthquake cycle time here is not 
known, although it is surely more than 160 yr. The fits are rather 
poor, as is perhaps attributable to the influence of the 
subparallel faults which essentially spread the plate boundary 
deformation over a much broader zone. Even if we focus only on the 
data points to the east of the San Andreas, the fits are still 
unsatisfactory. Clearly, a simple fault model with asymmetric 
displacement rates cannot adequately describe the surface 
deformation field in a region with active sub-parallel faults. 
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Conclusions 

We have presented a physical model which allows the prediction 
of surface strain rates and surface displacement rates at a 
strike-slip plate boundary. The major features of this model are 
that it involves relatively thin elastic lithospheric plates (17 km 
to 25 km), locked between earthquakes over shallow depths (9 to 11 
km) at the plate margin but with continuing aseismic slip below the 
locked zone, and that the surface plates are coupled by a Maxwell 
viscoelastic asthenosphere to upper mantle motions compatible with 
overall plate motions. The surface deformation is shown to be 
localized to the plate boundary (to within approximately three 
lithospheric plate thicknesses) due to deep aseismic slip below a 
locked brittle upper crust, modelled as stable shear sliding of 
crack faces at constant resistive stress in the anti-plane strain 
deformation of an edge cracked elastic strip. 

The viscoelastic coupling is accomplished approximately by a 
generalized Elsasser model which provides the time-dependence of 
plate loading as well as the time dependence of the surface 
displacement and strain rates following sudden slips of the locked 
zone. Thus the complete earthquake cycle is modelled as sudden 
slip on the shallow locked upper crust (which will partially 
penetrate the aseismic shear zone) during an earthquake, which 
loads the asthenosphere. This is followed by the relaxation of the 
asthenosphere which reloads the surface plate, and this loading 
causes continuing aseismic slip below the now locked zone. Because 
of the relaxation effect of the asthenosphere, the stress 
accumulation and the surface deformation rates are non-linear in 
time over the earthquake cycle. Finally it should be emphasized 
that while the San Andreas fault may accomodate an overall slip 
rate of 35mm/yr, the contemporary slip rate below the locked zone 
may be much smaller than 35mm/yr. The substantially higher slip 
rate shortly after an earthquake rupture makes up the deficiency 
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and produces an average rate of 35mm/yr over the earthquake cycle. 

The essential parameters in the model are the plate velocity 
(that part of which is accomodated by the San Andreas fault) , the 
locked zone depth, the surface plate thickness, the cycle time and 
the relaxation time of the modified Elsasser asthenospheric 
foundation. By arguments of previous geologic, geodetic and 

seismic studies, we have assumed the plate velocity to be 35mm/yr, 
the cycle time for the San Andreas as a whole to be 160yr, and the 
locked zone depth to be 9-llkm. Thus the only free parameters in 
the model are the lithospheric plate thickness and the relaxation 
time of the foundation. The strain rate data from various 

locations of the San Andreas fault collected by Thatcher (1983) 
provides the constraints to these parameters, and it is found that 
the plate thickness cannot deviate significantly from 17km to 25km, 
while the relaxation time for the foundation cannot deviate 
significantly from lOyr to 16yr. We conclude that the 
asthenosphere appropriate to describe crustal deformation on the 
earthquake cycle time scale lies in the lower crust and perhaps in 
the crust-mantle transition zone and has an effective viscosity 
between about 2xl0 18 and 10 19 Pa-s, depending on the thickness 
assigned to the asthenospheric layer. 

To gain further insight into the nature of the contemporary 
distribution of surface deformation at the San Andreas fault, the 
model with parameters constrained as described above has been used 
to predict the surface displacement rate and surface strain rate 
profiles at several locations along the San Andreas where such data 
has been obtained from repeated geodetic line measurements. These 
include the geodetic net from the Palmdale area, from the Point 
Reyes area, and from the Coachella Valley-Salton sea area. The 
results indicate that good fits to data can be obtained, although 
the simple single fault geometry as contained in the present 
version of the model clearly cannot account for data contaminated 
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by motions on branch faults or by material properties which are not 
distributed symmetrically about the San Andreas. 

The model is seen here to fit the principal types of geodetic 
data now available along .the San Andreas, consisting of strain rate 
near the fault trace as a function of time since the last great 
earthquake as assembled by Thatcher (1983), and of contemporary 
surface displacement rate as a function of distance from the fault 
trace. Thus, we conclude that such geodetic data does not require 
one to contemplate models involving an extremely thin elastic 
lithosphere that terminates jus't below seismogenic depths 
(Thatcher, 1983) . While Turcotte et al (1984) were able to fit the 
San Andreas data of fig. 6 for Y s as a function of time, using 
their four layer model, it is not yet clear whether their model 
also describes suitably the contemporary y s and u s as a function of 
distance from the trace. 

Even with its apparent successes, the model that we present 
here has evident limitations. For example, the aseismic shear zone 
at depth is assumed to creep-slip under constant resistive stress. 
However, in reality it may be expected that this zone supports 
larger stresses just after the earthquake in the shallow crust and 
then relaxes with time, hence producing a second source of time 
dependence which is not accounted for in the present model. It 
seems possible to incorporate this viscous shear effect of the deep 
aseismic shear zone approximately by modifying the boundary 
condition (5) from its present elastic relation between thickness 
average stress and slip to an appropriate viscoelastic relation. 
However, based on short term post-seismic deformation observation 
and modelling in subduct ion zones, Thatcher (1984) and Thatcher and 
Rundle (1984) suggested that the relaxation time for the deep shear 
zone is much shorter than for the asthenosphere, the former being 
on the order of a year or less and the latter tens of years . 
Hence, we may expect that our present results would not be affected 
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by this modification if one is mainly concerned with long term 
interseismic plate boundary deformation and geodetic data 
interpretation, as we are in this paper. Also, there are 
simplifying approximations inherent in our use of the Elsasser 
foundation concept and Turcotte-Spence edge-cracked strip solution 
which could be improved upon in future work. 
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Appendix A: Solution of differential equation 

We reproduce here the governing equation (4) 

(a + p3/9t) 9 2 u/9y 2 = 9u/9t (Al) 

for the part of u that must respond to the first periodic sawtooth 
term in (6) . This part is subject to the boundary condition (5) 

G 9u/9y = 2k (u - D/2) on y = O (A2) 

with D given by the sawtooth. Introducing the normalized 
parameters 

Y b y/7p , T m at/p = t/t r , 

X s G/2k-/p = ( 8 / 7E 2 ) In [1 / sin (7CL/2H) ] 

(Al) and (A2) become 

■ (1 + 9/3t) 9 2 u/9y 2 = 9u/9t (A3) 


and 

X 3u/3y - u = -D/2 on Y = 0 (A4) 

We first look for the solution to the fundamental problem of a 
period slip function D = D o e it0T and later superimpose all 
frequencies CO to form the complete solution to our problem. Let 
u(Y,T) = u o (Y)e it0T . Then (A3) and (A4) become 


(1 + iCO) 9 2 u 0 /9y 2 = iC0u o 


(A5) 
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X du o /dY - u o - - D o /2 on Y = 0 + (A6) 

from which it may be shown that 

u 0 (Y) = [D 0 /2B(C*»] exp (-M(CO)Y - i [N (CO) Y + \|/(C0) ] } 
where 

I 2 

M (£0) = VC0[p(C0) + (D]/2 [p (CO) ] 

N (CO) = ■>}<&[ p(C0) - CO] /2 [p (CO) ] 2 

p (co) = V : i + co 2 

B (CO) = yj 1 + 2XM(C0) + X, 2 CO/p (CO) 

\|/ ( CO) = arctan {A, N(CD)/[1 + X M(C0)]} 

Hence, for D(T) = Im (D Q e^ wt ) = D o sin COT, 

U (Y, T) - [D o /2B(C0)] exp [-M(0))Y] sin [COT - N(C0)Y - \|/ (CO) ] (A8) 

As the periodic sawtooth function may be represented by an infinite 
sine series, the solution for u(Y,T) also takes the form of a sine 
series, made up of fundamental solutions given by (A8) with the 
appropriate frequency CO. The complete form of u(Y,T) is given by 
equation (7) . For numerical convergence, we have removed terms 
summing to known discontinuous functions of time (e.g., sawtooths) 
from within the summation signs. The actual expression used for 
numerical computation is given below for the thickness averaged 
displacement rate: 
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1 5u , 1 e 

— ffft) = - + — 

2 l+A, 


v , at 

pi 


1 /I--M (_L + Y \ Is- I 

4 \ 2 T cy / \ 1 + X / t 2 


oo 

+ X (a n 8in(® n t/t r ) + |3 n COS (® n t/t r )] 
n-1 


(A9) 


where 


-M n Y 


a = 
n B 


e . e 

sm (vi/ + YN ) 

n n 1+Jt 


(i) (t^ + y ) 


-Y 


P = cos (VI/ + YN ) - — 

B n " n 1+X 


CD = 27tnt /T 

n r cy 


and where M n , N n , p n , B n , and V|/ n are the functions of CO given above 
evaluated at co = CD . 


Also, for thickness averaged shear strain rate, we use 


a 2 u 

- 1 ±L 

1 + ± 1 

l 

(— - y) 


dydt 

K 1+X, 

2 

\ K ) 

\ l+X / 

\ 2 T C J 


pi 


-Qtc 2 | d/6) - (t/T ) + (t/T ) 2 

cy cy 


[ 


] 


(A10) 


~ X [ 9 n Sin(a) n t/t r ) + \ C0S <®n t/t r ) ] 

n=l 
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where 



Notice that all the coefficients a n , P n , 0 n and y n of the sine and 
cosine terms in the summations vanish as 1/n 2 as n — » °°, which 
accelerates the series convergence. Equations (A9) and (A10) are 
the actual expressions used for numerical evaluation of equations 
(12) and (13) . 
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Appendix B: Accuracy of surface displa c ement estimation by eqn . 

(12) and (131 

We consider the deformation of an elastic plate of thickness H 
with an edge crack of length a (sH-L), loaded at the base of the 
plate by imposed displacement rate V pl /2 for y>0 and -V pl /2 for y<0, 
as shown in Fig. Bl. This loading configuration represents an 
extreme case of lithosphere/asthenosphere coupling where the 
elastic lithosphere is driven by a rigid foundation moving at 
uniform velocity consistent with the far field plate motion. We 
shall call this the uniform base driven model. 

For this extreme case, 2-D elasticity theory leads to a 
governing equation for the only non-vanishing displacement rate 
component w(y, z) : 


V 2 w = 0 (Bl) 

subject to the boundary conditions indicated in Fig. Bl. 

The solution to this boundary value problem has been obtained 
by means of the conformal mapping technique and the details are 
reported elsewhere (Li and Lim, unpublished manuscript, 1986) . 
Here we summarize the resulting expressions relevant to the present 
discussion : 

(1) Displacement rate in Body: 

p 

w(y, z) = (V /2) Re f . (B2) 

P j / 2 2 2 

O ) (1-k t ) 


where 


k 


sin a 
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Jt/2 


F = 


J 


d8 

I 2 2 ~ 

■y 1-k sin 0 


is the elliptical integral of the 1st kind 


p = Re 


y^ 


2 2 
tanh (y+iP) /tan a 


a = 7C a/2H; y ■ 7ty/2H; P = 7Cz/2H 


(2) Surface displacement rate (at z=H) 


u a (y) » w(y,z=H) = V pl /2 


1 - (1/F) 


A 

1 

1/k 


dt 




2 2 
1) (k t -1) 


(B3) 


WT7 


2 2 

where “k = -Jl + cot h y / tan a 


(3) Stress rate in body: 


X - i x = 

yx zx 


G(V pl /2F) (7C/2H) 


7 


sin[iY - (p-a) ] sin[iY - (P+a) ] 


(B4) 


(4) Surface strain rate (at z=H) 


(V /2F) (7C/2H) 

Y s (y) = 

2 2 
y smh Y + cos a 


(B5) 
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We compare predictions of surface deformation rate from the 

-E 

exact solution, here called u g (y), just obtained to that predicted 
from two versions of eq. (12) in the text: 

u*(y) = u(0 ) S (y) + [u(y) - u(0 )] (B6) 

and 

u®(y) = [0(0) /2k] S (y) + [u(y) - u(0 + )] (B7) 


where u(y) again denotes the thickness average of w(y,z) and o(y) 
is the thickness average of T„ v (y,z). These comparisons, calculated 

yx 

for L = 0.25H, 0 . 5H and 0.75H are shown in figure B2 . In equation 

(B6) and (B7), u and O are thickness averages of displacement and 

shear stress based on the exact solution and k is the stiffness of 

eq. (5) . When the edge cracked strip is loaded by surface 

tractions that act remotely from the crack section, 0(0) = 2ku(0 + ) 

so that versions A and B are equivalent. In the present case, 

however, owing to the imposition of displacements directly at the 

bottom of the crack, we find that the crack slip profile deviates 

to larger values at depth from the profile expected from remotely 

applied loading, so that u(o + ) slightly exceeds 0(0) /2k. Thus, 

0(0) is numerically very close to 2ku(0 + ) only at small L/H. The 

latter version B then more successfully embodies the notion that 

the deformation perturbation at the fault trace is governed by the 

net force transmitted across the locked ligament. We do not 

replace the second u(0 + ) in version B by 0(0) /2k, since that would 

give a discontinuity to u s (y) at y = 0, equal to the difference 

u (0 + ) - 0 ( 0 ) /2k. (The predicted u s (y) at y = 00 in version B is, 

however, too small by the same difference) . Note that the second 

• 

u(0) makes no contribution whatever to the surface strain rate y = 


\ 
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du s /dy, so that y in version B is determined entirely by the net 

• • 

force distribution, that is by the distribution a(y) = G du(y)/dy. 

The estimated surface strain rate profile may be obtained from 
a direct differentiation of (B7) with respect to y: 

Y s (y) = [C(0)/2k] dS(y)/dy +du(y)/dy (B8) 

This may be compared with the exact solution (B5) and the 
comparison is shown in Figure B3, for three different locked zone 
depths. Furthermore the ratio of the surface strain rate at the 
fault trace for the exact solution to that for the approximate 
solution B is shown in Fig. B4, as a function of locked zone depth. 

The surface displacement and strain rates for the edge-cracked 
strip can be regarded as the response of an elastic body to the 
base loading rate T x . The heavy solid curve in Figure B5 shows the 
distribution of loading rate for the base driven model discussed in 
this appendix, when L/H = 0.5. The loading rate is concentrated 
near the fault zone and thus introduces the moderate errors shown 
in figures B2 to B4 . The loading rate distributions at various 
fractions of the cycle time are also shown in figure B5 for the 
model as developed in the text, using t r /T cy = 12 yr/160yr = 0.075. 
Evidently, those loading rates has, by comparison, a very broad 
distribution at mid- and late-cycle times creating base loading 
conditions that are more compatible with the approximations in our 
modelling procedures. 
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Table 1: Calculated and Observed Strain Rate in the Palmdale Area 



y 

(Ay) rms 

Y (a) 

Y (b) 

Y(c) 

Y (King & Savage 
1984) 

+Palmdale 

-0.3 

3.6 

0.42 

0.41 

0.38 

0.37 ± 0.02 

(1971-82) 

+S.A. Region 

4.1 

10.4 

0.37 

0.38 

0.33 

0.34 ± 0.01 

(1973-83) 

+Tehachapi 

17.0 

20.2 

0.29 

0.31 

0.25 

0.21 ± 0.01 

(1973-83) 

Los Padres 

-19.7 

19.7 

0.27 

0.29 

0.23 

0.21 ± 0.02 
(McGarr et al., 
1982) 

+Sub Garlock 

31.0 

21.1 

0.21 

0.23 

0.18 

0.17 ± 0.02 


(1973-83; called 
Garlock region) 


+ subdivision of Tehachapi net 

Y in (jistrain/yr) and y, (Ay) rmg in km, positive to the NE. All calculations made 
for t = 123 yr, T cy = 160 yr and (a) V pl = 35 mm/yr, L - 9 km, H = 22.5 km, t r = 12 
yr; (b) V . = 35 mm/yr, L = 11 km, H = 20 km, t r = 14 yr; (c) V . = 32 mm/yr, L = 9 
km, H = 25 km, t r = 12 yr. 
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Figure 1 : 


Figure 2 : 


Figure 3 : 


Figure Captions 

(a) Elastic lithosphere coupled to a viscoelastic 
asthenosphere driven by deep mantle movement. The 
shaded area indicates the shear zone sliding at 
constant resistive shear stress below the locked 
brittle zone, (b) A cross-sectional view of the 
lithosphere at the plate boundary, modelled as an 
edge cracked strip in anti-plane strain. 

Step movements at plate boundary decomposed into 
steady motion and a time periodic sawtooth motion, 
with cycle time T . 

Surface shear strain rate y s / (V pl /H) as a function of 
distance from the fault, at various fractions (0.0, 
0.1, 0.2, 0.3, 0.4, 0.6, 0.8, 1.0) of an earthquake 
cycle T cy . These figures are based on T cy = 160 yr., 
t r = 12 yr, and with contrasting locked depths of (a) 
L = H and (b) L = 0.3H. Note that in case (b) , 
where shear slip occurs at depth, the shear strain 
rate is much more localized at the plate boundary, 
as is particularly evident during later parts of the 
earthquake cycle. 

Displacement of shear crack wall below locked zone 
as a function of depth, at various fractions of 
cycle time. Same parameters as for Figure 3b. 


Figure 4 : 
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Figure 


Figure 


Figure 


Figure 


Figure 


Figure 


5: Map of California with locations mentioned in this 

paper. SC - Shelter Cove; PA - Pt . Arena; FR - Ft. 
Ross; PR - Pt . Reyes; SF - San Francisco; CP - 
Carrizo Plain; MD - Mojave Desert; PD - Palmdale; CV 
- Coachella Valley. (After Allen, 1968; Thatcher, 
1983) . 

6: Shear strain rate decay with time in an earthquake 

cycle. The data points with error bars are from 
Thatcher (1983). The solid circles are from 
northern California along the 1906 rupture zone of 
the San Andreas fault. The open circles are from 
Southern California along the 1857 rupture zone. 
Curve fits are based on model with indicated 
parametric values. 

7: Displacement rate profiles, calculated for the 

parameters indicated. 

8: Geodetic networks in the Palmdale area (after. King 

and Savage, 1984) . 

9: Comparisons of theoretically predicted displacement 

rate profiles to geodetic data from King and Savage 
(1984). Stations, in order of increasing y, are: 
Mt . Pinos, Frazier, Sawmill, Tejon 41, Tecuya, 
Police, Thumb, Wheeler 2, Diorite, Gneiss and Tejon 
32 . 

10: Geodetic networks in the Point Reyes - Santa Rosa - 

Napa area (after Prescott and Yu, 1986) . 
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Figure 11: Comparisons of theorectically predicted displacement 

rate profiles to geodetic data from Prescott and Yu 
(1986), with error ranges indicated. Circle symbols 
are data from Pt . Reyes, Santa Rosa and Napa 
networks . Triangle symbols are from Geyser network 
further north. The two data points at approximately 
-20 km and -40 km are associated with the Pt. Reyes 
Head and Farallon Islands markers. 

Figure 12: Geodetic networks in the Coachella Valley - Salton 

Sea area (after Savage, 1983) . 

Figure 13: Comparisons of theoretically predicted displacement 

rate profiles to geodetic data. Data from Savage 
(1983) . Circle symbols are data to the NW of the 
San Andreas fault . Triangle symbols are data to the 
SW of San Andreas fault . 

Figure B1 Base driven model of lithosphere plates loaded by 

uniform velocity V pl /2 at y>0 and -V pl /2 at y<0. 

Figure B2 Comparison of surface displacement rate profiles for 

base driven model between exact (eqn. B3) and 
approximate (eqn. B6 and B7) solutions, for (a) L = 
0.25H (b) L = 0.5 H and (c) L = 0.75H. 

Comparison of surface strain rate profiles for based 
driven model between exact (eqn. B5) and approximate 
(eqn. B8) solutions, for (a) L = 0.25 H, (b) L= 0.5 
H and (c) L = 0.75 H. 


Figure B3 
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Figure B4 Ratio of surface strain rate at fault trace between 

exact and approximate solution, as a function of 
locked zone depth. 

Figure B5 Loading rate distributions on base of edge cracked 

strip, for base driven model of figure Bl and for 
viscoelast ically coupled model developed in main 
text . 
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